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Abstract 

In previous work by Y. Guo and G. Rein, nonlinear stability of 
equilibria in stellar dynamics, i.e., of steady states of the Vlasov- 
Poisson system, was accessed by variational techniques. Here we pro- 
pose a different, non- variational technique and use it to prove nonlin- 
ear stability of the King model against a class of spherically symmetric, 
dynamically accessible perturbations. This model is very important 
in astrophysics and was out of reach of the previous techniques. 



1 Introduction 

In astrophysics a galaxy or a globular cluster is often modeled as a large 
ensemble of particles, i.e., stars, which interact only by the gravitational field 
which they create collectively, collisions among the stars being sufficiently 
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rare to be neglected. The time evolution of the distribution function / = 
f{t,x,v) >0 of the stars in phase space is then given by the Vlasov-Poisson 
system: 

dtf + v-VJ-V.U-\/J = 0, (1.1) 



AU = Anp, lim U{t,x) = 0, (1.2) 



p{t,x)= J f{t,x,v)dv. (1.3) 

Here tGM denotes time, and x,fGM^ denote position and velocity, U = 
U{t,x) is the gravitational potential of the ensemble, and p = p{t,x) is its 
spatial density. 

A well-known approach to obtaining steady states of the Vlasov-Poisson 
system is to make an ansatz of the form 

fo{x,v) = (P{E), 

where the particle energy E is defined as 

E:=^\v\^ + Uo{x). 

Since for a time-independent potential Uo = Uo{x) the particle energy is con- 
served along the characteristics of the Vlasov equation (jl.lll it remains to 
make sure that for the chosen the resulting semilinear Poisson equation 

A?7o = 47r j 0Q|t;|2 + f/o^ dv (1.4) 

has a solution. For a very large class of ansatz functions this approach leads 
to steady states with finite mass and compact support, cf. |RRj . Steady states 
where as above the distribution of the particles in phase space depends only 
on the particle energy are usually called isotropic and are always spherically 
symmetric. A central questions in stellar dynamics, which has attracted 
considerable attention in the astrophysics literature, cf. |BT| IFPj and the 
references there, is the dynamical stability of such steady states. 
The Vlasov-Poisson system is conservative, i.e., the total energy 

nif):=EUf)+Epotif) = \ j \v\^f{x,v)dvdx~^ j \VUf{x)\^dx 
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of a state / is conserved along solutions and hence is a natural candidate for 
a Lyapunov function in a stability analysis; Uf denotes the potential induced 
by /. However, the energy does not have critical points, i.e., the linear part 
in an expansion about any state /o with potential Uq does not vanish: 

H(/)=H(/o) + jj (^l^l' + f^o) if-fo)dvdx-^ j \WUf-VUo\''dx. 

To remedy this situation one observes that for any reasonable function $ the 
so-called Casimir functional 



C{f):= J J ^f{x,v))dvdx 
is conserved as well. If the energy-Casimir functional 

nc:=n+c 

is expanded about an isotropic steady state, then 
ndf) = nc{fo) + j I {E + <^>'{fo)){f-fo)dvdx 

/ mf-VUo\'dx + ^ J I <^"{h){f-hfdvdx+... 
=: Hcih) + DHcih) [f - h] + D^Hcih) [/ - /o] + • • • • 

Now one can try to choose $ is such a way that at least formally /q is a 
critical point of the energy-Casimir functional, i.e., $'(/o) = —E. If 0' < on 
the support of /o the former relation holds on this support if $' = —0^^. The 
formal second order variation in the above expansion then takes the form 

D'nc{h)\g] = \ [ [ —^ydvdx-^ [ \VMg\'dx. (1.5) 

It is natural to expect that positive definiteness of this quadratic form should 
imply stability for /g. Ever since the seminal work of Antonov |Anj there 
have been vigorous efforts in the astrophysics community to establish this 
positive definiteness and to derive stability results in this way. 

An important step in this direction was to show that the above quadratic 
form is positive definite on linearized, dynamically accessible perturbations. 
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To make this precise we define the Lie-Poisson bracket of two functions /i,/2 
of X and v as 

{/i,/2}:=V./i-V./2-V./i-V./2. (1.6) 
Then the following holds: 

Lemma 1.1 Let 0'<O and let /iGC^(M^) be spherically symmetric with 
support in the set {/o>0} and such that h{x,—v) = —h{x,v). Then 

D'ncm{foM]>-\ [ <P\E) 

^ Jfo>0 

Here Uq denotes the radial derivative of the steady state potential. Since 
Uq is radially increasing, the right hand side in the estimate above is indeed 
positive for hy^O. We refer to |KSt ISDLPj for astrophysical investigations 
where this result is used to analyze linearized stability. We do not go into 
the reasons why perturbations of the form {fo,h} are called dynamically 
accessible for the linearized system, but for the sake of completeness we 
provide a proof of this elegant result in the Appendix. Despite its significance 
the result is still quite a distance away from a true, nonlinear stability result. 
There are at least two serious mathematical difficulties. Firstly, it is very 
challenging to use the positivity of D'^T-Ccifo)[9] to control the higher order 
remainder in the expansion of the energy-Casimir functional |Waj . This is due 
to the non-smooth nature of /o = 4>{E) in all important examples. Secondly, 
even if one succeeds in controlling the higher order terms, the positivity of 
D'^T-lc{fo)[g] in the lemma is only valid for certain perturbation of the form 
g = {fo,h}. This class of perturbations is invariant under solutions of the 
linearized Vlasov-Poisson system, but it is not invariant under solutions of 
the nonlinear system. 

To overcome these difficulties a variational approach was initiated by 
WOLANSKY |Wolj and then developed systematically by Guo and Rein 
[HTl IHrTI inR2l IHMl IHrI eh Ea EG] . Their method entirely avoids 
the delicate analysis of the second order term D'^Hcifo) in (jl.5|) . and it 
has led to the first rigorous nonlinear stability proofs for a large class of 
steady states. More precisely, a large class of steady states is obtained as 
minimizers of energy-Casimir functionals under a mass constraint / / = M, 
and their minimizing property then entails their stability. In particular, all 
polytropes fo{x,v) = {Eq — E)'^^ with 0<fc<7/2 are covered; here -E'o<0 is 
a certain cut-off energy, and (■)_). denotes the positive part. For k>7/2 the 



\x-v\ 



E, 



x-v 



+ -U',h' 
r 



dvdx. 
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corresponding steady state has infinite mass and is therefore unphysical. In 
addition, many new stable galaxy models were established. The variational 
method has also been investigated in iHl iLMRl 1^ IWo2] . 

Despite its considerable success, the variational approach has drawbacks 
and limitations, the main one being that by its very nature it can not access 
the stability of steady states which are only local, but not global minimizers 
of the energy- Casimir functional. Since the existence of the steady state as 
a (global) minimizer is aimed for, certain growth conditions on the Casimir 
function $ are needed, which are not satisfied for all steady states with 0' < 0. 
Most notably, the King model obtained by 



is the single most important model which is currently out of the reach. It 
describes isothermal galaxies and is widely used in astrophysics. The corre- 
sponding Casimir function (|2.2j) has very slow growth for f—*oo, and as a 
result the variational method fails. 

The aim of the present paper is to develop a new approach to nonlinear 
stability results for steady states which need not be global minimizers of the 
corresponding energy-Casimir functional by exploiting Lemma ITTTl Although 
we are aiming for a general approach, we focus here on the King model and 
as a first step establish its nonlinear stability against spherically symmetric, 
dynamically accessible perturbations. 

The paper proceeds as follows. In the next section we formulate our re- 
sults. The nonlinear stability of the King model is an easy corollary of the 
following main theorem: In a certain neighborhood of the King model the 
potential energy distance of a perturbation can be controlled in terms of 
the energy-Casimir distance. In particular, within a certain class of pertur- 
bations, which is invariant under solutions of the nonlinear Vlasov-Poisson 
system, the King model is a local minimizer of the corresponding energy- 
Casimir functional. The resulting stability estimate is more explicit that the 
ones obtained by the variational approach. The main part of the work is then 
done in Section 3 where the local minimizing property of the King model is 
established. In an appendix we give a proof of Lemma 11.11 
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2 Main results 



We start with a steady state /o with induced potential Uq and spatial density 
Po, satisfying the relation 

fo{x,v) = ME):={e'''-''-l)^, E ■.= ^\v\' + Uo{x). (2.1) 

The cut-off energy ii^o<0 is a given negative constant and (■)+ denotes the 
positive part. The corresponding Casimir function in the sense of the intro- 
duction is 

*o(/):=(l + /)ln(l + /)-/. (2.2) 

The existence of such King models, i.e., of suitable solutions of the resulting 
semilinear Poisson equation ()1.4p . is established in jRRj . Such a model has 
compact support 

supp/o = {{v,v) G I E{x,v) < Eo} =: {E < Eo}, 

and it is spherically symmetric. A state / is called spherically symmetric if 
for any rotation A G SO (3), 

f{x,v) = f{Ax,Av), x,veR^. 

It is well known that non-negative, smooth, and compactly supported ini- 
tial data /(O)gC^(M^) launch unique global smooth solutions t^f{t) of 
the Vlasov- Poisson system |PH ILP| ISchj . If the initial datum is spherically 
symmetric then this symmetry is preserved, and the modulus of the particle 
angular momentum squared, 

T I |2 Ii2ii2 / \2 

L:=\xxv\ =\x\ \v\ —[x-v) , 

is conserved along characteristics of the Vlasov equation. Hence for any 
smooth function $ such that $(0,L) = 0, L>0, the functional 

j I ^{f,L)dvdx 

is conserved along spherically symmetric solutions of the Vlasov-Poisson sys- 
tem; unless explicitly stated otherwise integrals J always extend over R^. 
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We consider the following class of perturbations: 
Sfg := | / G L^(R^) I / > spherically symmetric, 

$(/,L)= [ I <l>(/o,L) for all $€^^([0,00^ with 



$(0,L) = 9/<I>(0,L) = 0, L>0, axidd)^ bounded}. 



As noted above, the class 4S/QnC^(]R^) is invariant under solutions of the 
Vlasov-Poisson system. Moreover, functions in Sf^ are equi- measurable to 
/o, i.e., for every r > the sets {/ > r} and {/o > r} have the same measure, 
in particular, ||/||p= ||/o||p for any L^'-norm, pG [l,C)o]. 

For the Casimir function $o defined in (12. 2j) we define the energy- Casimir 
functional as in the introduction. Then for /GiS/q we have 



dvdx 



nc{f)-nc{fo) = jj [%{f)-%ifo) + iE-Eo)if-fo)] 

-^j \^Uf-VUo\^dx; 

notice that J f = f fo which allows us to bring in the term Eo{f — fo). Now 

(E-Eo)(/-/o)>-$'o(/o)(/-/o) 
with equality on the support of fo, and hence for /g5/q, 

Mf)-Mfo) + {E-Eo){f-fo)>l^^ inf „ <^>o(r)(/-/of 

^ 0<T<||/o||cx) 

>Co{f-foy (2.3) 

where Cq := l/(2 + 2| |/o| U); notice again that 1 1/| U = I |/o| U for any feSf^. 
The deviation from the steady state is going to be measured by the quantity 

d{f,fo) ■■= jj [Mf)-Mfo) + {E-Eo){f-fo)]dvdx 

+ ^ [ \VUf-VUo\^dx, (2.4) 



57r „ 

which, as we have seen, controls both ||/ — /0II2 and ||Vt// — Vf/o||2, and 
satisfies the following relation to the energy-Casimir functional: 

d{f,fo)=nc{f)-nc{fo) + ^ j \vUf-vUo\'dx. (2.5) 

Our stability result is the following: 
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Theorem 2.1 There exist constants S>0 and C>0 such that for any solu- 
tion tt— of the Vlasov-Poisson system with /(O) g5/„ nC^(M^) and 

d{f{0)Jo)<6 

the estimate 

difit)Jo)<CdifiO),fo) 

holds for all time t>0. 

Remark: In order to better understand the perturbation class 5/^ we show 
that it contains spherically symmetric, dynamically accessible perturbations 
by which we mean the following: Let H = H(x,v) &C'^{M.^) be spherically 
symmetric, and let g = g{s,x,v) denote the solution of the linear problem 

dsg + V„H-V^g-V,H-V,g = 0, i.e., dsg{s) = {H,g{s)}, 

with initial datum g{0) = fo', we assume that H is such that this solution 
exists on some interval / about s = 0. Then for any sG/, g{s) is spherically 
symmetric and equi-measurable with /q. Moreover, for any function $ as 
considered in the definition of 5/^, 

ds^gis),L)=df^gis),L){H,gis)}, 
and by a simple computation, {H,L} = 0. Hence 

j ^{g{s),L)dvdx = J I {df^{H,g{s)} + dL^{H,L})dvdx 

= jJ{H,^{g{s),L)}dvdx = 

after an integration by parts. This means that g{s) GiS/g for any such gen- 
erating function H and any s. The only undesirable restriction in the class 
Sfg or in the generating functions H respectively is the spherical symmetry 
which hopefully can be removed in the future. 

The stability result Theorem 12.11 is easily deduced from the following 
theorem: 

Theorem 2.2 There exist constants So>0, and Cq > such that for all f G 
iS/q with d{f,fo) <6o the following estimate holds: 

nc{f)-ncUo)>Co\\vUf-vUo\\l 
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Before going into the proof of this theorem, which will occupy the rest of this 
paper, we conclude this section by deducing our stability result from it. 
Proof of Theorem mU Let 5 ■. = 5o{l + l/ {AtiCq))-^ with So and Co from 
Theorem 12 .21 Consider a solution f{t) of the Vlasov-Poisson system with 
f{Q)eSf,nCl{R^) and 

d{f{0),fo)<S<6o. 
Then by continuity we can choose some maximal t* g]0,cxd] such that 

d(/(t),/o)<5o, te[o,r[. 

Now /(t)GiS/o for all t, and hence Theorem 12.21 the relation ()2.5|) of d to 
the energy- Casimir functional, and the fact that the latter is a conserved 
quantity yield the following chain of estimates for tE [0,t*[: 

d{f{t)jo)=ncif{t))-nc{fo) + l^\mj^t)-'^Uo\\l 

< nc{f{t))-nc{fo) + ^ {nc{f{t))-nc{fo)) 



< ( 1 + H(/(0),/o)<5o. 



This implies that t* = oo, and Theorem 12. II is established. □ 



3 Proof of Theorem 

Theorem 12.21 is proven by contradiction. There are two main ingredients: 
The first part (Subsection 13. Ij) is a general argument to establish that if 
the estimate in the theorem fails, then there exists a non-zero function 
g such that D'^'Hc{fo)[g]<0- The second (Subsection 13. 2j) is to use the 
measure-preserving property incorporated in our perturbation class 5/„ to 
conclude that g = {fo,h} for some function h. This leads to a contradiction 
to Lemma fTTT] f Subsection . 

3.1 Existence oi gj^O with D^nc{fo)[g]<0 
The aim of this subsection is to prove the following result: 
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Lemma 3.1 Assume that Theorem \2.^ were false. Then there is a function 
(7G-L^(]R^) which is spherically symmetric, supported in the set {E<Eq}, 
even in v, i.e., g{x,—v)=g{x,v), and such that 

^\mx=h (3.1) 



D'nc{fo)[g] = l [ [ ^o{fo)9'dvdx-l<0, (3.2) 

^ J J{fo>0} 

J j df^fo,L)gdvdx = (3.3) 



for all functions $ as specified in the definition ofSfg. 

Proof. If Theorem 12.21 were false, then for any riGN there exists fnESf^ 
such that 

d{fnJo)<- 
n 

but 

nc{fn)-nc{fo) < ^||Vf//„ - V?7o||i (3.4) 
in particular, fny^fo- We define 

a„:=^||Vf//„-Vf/o||2, gn:= — {fn-fo) (3.5) 



SO that 



^"Vf/,J|^ = 1 (3.6) 



»7r 

and fn = fo + CTngn- By (jTH) and dSH)), 

a^<c^(/n,/o)<^. (3.7) 
Bounds on (gn) and weak limit. 

As a first step in the proof of the lemma we need to establish bounds on the 
sequence ((?„) which allow us extract a subsequence that converges to some 
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g which will be our candidate for the function asserted in the lemma. By 
dSini), dZll), dZni), and (jS3D we find that 

■\ I jl'^oifo + (Tngn) - $o(/o) + {E - Eo)fT„^„] dvdx-1 

= ^(^c(/.)-Hc(/o))<^^||Vt/,„-Vt/o||^ = ^. (3.8) 

Recalling the estimate (j2.3|) which is applicable since fn&Sf^ this implies 
that 

^ + ~;>~2 j I [■■■]dvdx>Co^ j I \fn- fo\'^dvdx = Co j j \gn\^dvdx, 

which means that the sequence is bounded in L^(]R^). Moreover, since 
the integrand [. . .] in ()3.8|) is non- negative we find that 

IH — > / / [. . ]dvdx> — I I {E — EQ)gndvdx, 

n ^nJ J{E>Eo} ^nJ J{E>Eo} 

where we used the fact that on the set {E > Eq} the steady state distribution 
/o and hence also $o(/o) vanish while $o(/n) >0. By ()3.7|) . 

/ / {E-Eo)gndvdx<2an-'0, n^oo. (3.9) 

J J{E>Eo} 

Now fix any Eq< Ei<0. Since lim|^|_,oof^o(3^) = it follows that E = 
E{x,v) > El for X or f large so that the set {E < Ei} gM.^ is compact. Hence 
(gn) is bounded in L^{{E<Ei}). In addition 

f f E-Eo 
gndvdx< i I — — gndvdx ^0; (3.10) 

{E>Ei} J J{E>Eo} El — Eq 

notice that (^n > outside the support of /o, i.e., on the set {E > Eq}. Thus we 
have shown that (gn) is bounded in L^nL^(M^). We extract a subsequence, 
denoted again by (gn) such that 



Qn^ 9 weakly in L 
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Since gn>0 on {E> Eq} and since (|3.1(J|) holds for any Eq<Ei<0 we con- 
clude that g vanishes a. e. outside the set {E<Eq} as desired. Since the 
functions Qn are spherically symmetric so is g. 

Proof of / TO) . 

In order to pass the week convergence into (j3.6|) we need better bounds for 
the sequence (gn)- Indeed, we can bound its kinetic energy: 

-\v\'^ + Uo{x) — Eq\ gndvdx = / / {E — Eq) gndvdx ^0 



{E>Eo} \2 / J J{E>Eo} 



by dSSI) so that 

-\v\'^ gn dvdx 

^{E>Eo} 2 



// 

J J {I 



is bounded, while 

/ / \\v\^ gndvdx<{EQ-UQ{Q)) 1 1 gndvdx 

J J{E<Eo} ^ J J{E<Eo} 

is bounded as well; recall that Uq is spherically symmetric and radially in- 
creasing. Now well known interpolation arguments imply that the sequence 
of induced spatial densities {pg„) is bounded in nL''/^(R^), cf. |R3t Ch. 1, 
Lemma 5.1], so without loss of generality, 

Pg^^Pg weakly in ^/^(M^). 

Let again Eq< Ei<0 be arbitrary but fixed and choose i?i > such that 
Uo{Ri) = Ei which implies that E{x,v) >Ei for |x| >-Ri. Then by (|3.1U|) . 



\PgJdx< / / gndvdx ^0. 

{\x\>Ri} J J{E>Ei} 



The fact that the sequence {pg„) remains concentrated in this way gives the 
desired compactness: 



VUg^^VUg Strongly in ^^(R^), 



cf. |R,3[ Ch. 2, Lemma 3.2]. Hence we can pass to the limit in ()3.(H) and find 
that g satisfies the condition p.ip in the lemma. 



12 



Proof of arm . 

Since \\gn\\2 is bounded it follows from (j3.7|) that 

\\(yngn\\2<Can^Q- 

Therefore, after extracting again a subsequence, (Tngn almost everywhere. 
By Egorov's Theorem there exists for every m G N a measurable subset Km C 
{E<Eq} with the property that 

vol {{E < Eq} \ Km) < — and lim OnQn = uniformly on Km', 

note that the set {E<Eq} has finite measure. In addition we can assume 
that KmCKm+i, meN. On the set {E<Eq}, 

[Mfn)-Mfo) + {E-Eo)angn] = l^'o{fo)Kgnf + ]:<^'^{fo + ran^^^^ 

for some r G [0, 1]. Since both /o and /o + 0"n5'n = fn are non- negative the same 
is true for fo + TO'ngn, and we can use the estimate 

IW)l = (Y^<i, />o 

the estimate ()3.8|1 . and the fact that the integrand [. . .] in ()3.8j] is non- negative 
by ()2.3|) to conclude that 



// \%{k)\9n?dvdx = \ f f [...]dvdx 



2 / / g*o'(/o + 'rcr„5(„)((T„5(„) dvdx 

< l + - + sup|(T„5f„| I I \gn\'^dvdx. 



Taking the limit n — >■ oo in this estimate we find that 



// lmfo)g'dvdx<i. 



If we observe the choice of the sets K^^ let m^oo, and recall the fact that 
(7 = outside the set {E<E^ the proof of (j3.2|) is complete. 
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Proof of / TO) . 

To prove (|3.3|) . the measure-preserving property of the set Sf^ plays the 
crucial role. Let $ = $(/,L) be a function as specified in the definition of 
that set respectively in the lemma. By Taylor expansion with respect to the 
first argument, 

<^{fn,L)-<^{fo,L)=df<!>{fo,L)angn + ^d}<!>{fo + Tangn,L){(Tngn)'^ 
for some rG [0,1]. If we integrate this identity and observe that, since G 

j I ^{fn,L)dvdx = j I <!>{fo,L)dvdx, 

it follows that 

df^{fo,L)gndvdx = -^an j ^ d^f^{fo + rangn,L) g^^dvdx -^0. 

As to the latter limit we note that is bounded, [gn) is bounded in L^(M^), 
and >0. On the other hand df^{fo,L) is supported on the compact set 
{E<Eo} and hence bounded. Since gn^Q in L^(M^), the identity ()3.3p 
follows as n — > oo. 

Conclusion of the proof of Lemma \S. 11 

The function g constructed above has all the properties required in the 
lemma, except that it need not be even in v. Hence we decompose it into its 
even and odd parts with respect to v: g = geven + dodd- We claim that ()3.1|) . 
fl3.2p . and ()3.3|) remain valid for the even part. Since 



Pg = J gdv = J C/evenC?t^ = Pgevcn 

we have VUg^^^^ = VUg, and ()3.1|1 remains valid. Since $o(/o) is even in v, 

1 > J l^'^oifo){gcven + goddYdvdx = J |"^$o(/o) {(gcvcn)^ + igoddf) dv dx 
- jj^^o(fo)ige^enfdvdx, 

i.e., (|3.2p remains valid. Finally, let $ be as in the definition of the set Sf^. 
Then df^{fo,L) is even in v so that in (j3.3p the odd part of g drops out, and 
the assertions of Lemma f3. II hold for gevcn- n 
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3.2 Characteristics and p = {/o,/i} 



In this subsection we construct a spherically symmetric function h such that 
g — {fo,h}. To this end we need to introduce variables which are adapted to 
the spherical symmetry: 



x-v 



r := \x\, w := 



, L:=\x X 



w is the radial velocity, and L has already been used above. Any spherically 
symmetric function of x and v such as the desired h can be written in terms 
of these variables, so h = h{r,w,L). Then 



{fo,h} = <l>oiE){E,h} = -<l>oiE) 



wdr+{^-U;,{r))d, 



and the equation we wish to solve for h reads 



--C/^(r) a, 



h- 



(3.11) 



In order to analyze its characteristic system 

. L 



r — w, w- 



we define for fixed L > the effective potential 

L 



^L(r):^Uo(r) 



2r2 



and observe that in terms of the variables r,w,L the conserved particle energy 
takes the form 

E = E{x,v) = E{r,w, L) = ]^w'' + *i(r) ; 

L only plays the role of a parameter here since L = 0. We need to analyze 
the effective potential The boundary condition for Uq at infinity implies 
that limj.^oo^'x,(r) = 0, and clearly, \im.T.^Q^ L{r) = oo. Moreover, 

U',{t) = % rp,{s)s'ds=:'^>0, r>0, 
r Jo ^ 
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since Uq is increasing, thus f/o(0) <Eq, and by (j2.ip . po(0) >0. Now 

^' (r) = ^ =0 ^ mo(r)-- = 0, 

and since the left hand side of the latter equation is strictly increasing for L > 
with limr^oo('^o('") ~-^/''') =^> and limr^olmQ^r) — L/r) = —(x there 
exists a unique > such that 

\&^(ri) = 0, \E'i,(r) < for r < r^, \E''i(r) > for r > r^. 

Moreover, since 



^^p(^) - = 47rrVo(r) + ^ > 

the implicit function theorem implies that the mapping ]0, oo[3 Lh^r^ is 
continuously differentiable. For r = rL, 

^'/(r^) = -2^^ + 3^ + 47rpo(r) = 47rpo(r) + ^ > 0. (3. 12) 

The behavior of \E'l implies that for any L > and \1'l(^l) < < there exist 
two unique radii 0<r_{E, L) <ri<r^{E,L) <oo such that 

^L(r±(E,L)) = E, and *L(r)<E ^ r_{E,L) <r <r+{E,L). 

Since \l/'^(r)7^0 for r^ri it follows again by the implicit function theo- 
rem that the mapping {E,L)\-^r±{E,L) is on the set {{E,L) eMx]0,oo[| 
^i(ri)<E<0}. 

Let r 1-^ (r(r),w(r),L) be a characteristic curve with L>0 and running 
in the set {E < Eq}, so that ^L(ri) < E < < 0. Then 



r(r)G[r_(E,L),r+(E,L)], u;(r) = ±v/2E-2^i(r(r)). 
By the equation (j3.11|) which we want to solve, 

-^/i(r(r),w(r),L) = -—^g{r{T),w{T),L), 

which we rewrite in terms of the parameter r as 
d 1 



—h{r,wir),L) = -————g{r,w{r),L), w;(r) = ±^/2E^^2¥^. 
dr (pQ[E)w{r) 
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Hence for {r,w,L) G {E < Eq} with L > we define h as follows: We let E :- 
^w'^ + '^Lir) so that r_{E,L) <r <r+{E,L), and 



<t>'o{E) J,_(E.L) V2S-2*i(.s) ^ 

Outside the set {E<Eo} we let /;, = 0. We need to make sure that we can 
consistently set h{r,0,L) = 0, i.e., the integral in the above definition must 
vanish for r = r^{E,L). To this end let $ be as specified in the definition of 
the class Sf^^. We want to apply the change of variables {x,v) i-^ {r,w,L) 
{r,E,L) to the integral in ()3.3|) . Now 

2 7 7 7 2 drdEdL 

dxdv = 8n drdwdL = 8n — , (3.14) 

v/2E-2M/^(r)' 

where we note that g is even in v and hence in w so that we can restrict the 
integral to w>0. We obtain the identity 

/» /• /"OO /"OO roo 

0= df^{fo,L)gdvdx = 8TT^ / / df^fo,L)gdrdwdL 

J JMJr-(E,L) A/2£;-2Wi(r) 



where M := {(i?,L)(x,f )|/o(x,f ) > 0}. The class of test functions 
df^{(f)Q{E),L) is sufficiently large to conclude that for almost all E and L, 



^+(^'^) g{r,^2E-2^L{r),L) ^^ ^ ^ 
r.(E,L) V2E-2^L(r) 



as desired. 



3.3 Contradiction to Lemma 11.11 

As defined above, h need not be smooth or even integrable, so in order to 
derive a contradiction to Lemma fl. II we need to regularize it. In order to do 
so it should first be noted that h as defined in ()3.13|) is measurable, which 
follows by Fubini's Theorem and the fact that by the change of variables 
formula the function 



g{s,^2E-2^L{s),L) ^ 

'^'^)^ /^j^ ^^j, I f l[*i(ri),£o]l^)l[r--(i?,L),r+(i?,L)](s)l[0,r](s) 
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is integrable; r <max{|x| | {x,v) Gsupp/o}. 
The cut-off h. 

As a first step in regularizing h we define for m large the set 

nm-=\{x,v)eR^\E<Eo-—, L> — 
I mm 

We want to approximate h by HIq^ . In order to analyze this approximation 
the following lemma will be useful: 

Lemma 3.2 There exists a constant Cm>0 such that for L>— and 
^LirL)<E<Eo, 

< Cm - 



^iE,L) ^/2E-2^l{ 

Proof. We first establish the following auxiliary estimate: For all m G N 
there exists a constant rim>0 such that for all L>l/m, "^l^tl) < E <Eq, 
andrE[r.{E,L),r+{E,L)], 

>Vm. (3.15) 



To see this let mGN and L,E,r be as specified. Then 

Eo>E>^L{r) = Uo{r) + ^>Uo{0)+ ^ 



2 



2r^ 2mr 

and hence r>{2m{Eo — Uo{0))^^^'^. Let i?: = max{|x| | (a;,f ) G supp/o}, i.e., 
Uq{R) = Eq. Then 

L < 2r^ {Eo - Uo{r)) < 2R' {E, - f/o(0)) . 

Hence if we assume that fl3.15|) were false, we can find a sequence {rn,Ln) — *• 
(f,L) in the set [{2m{Eo-UoiO)y^/^,R] x [l/m,2R'^ {Eo-Uo{0))] such that 

V*i..(''.,))-4'i,(riJ 



If r 7^ it follows that \E'i(r) > "^liri) and ^E'j^(r) = which is a contradiction 
to the uniqueness of the minimizer of "^i- So assume that f = ri. Now 



18 



recall that "^'j^ (fLn) = 0. By Taylor expansion at r = r^^ we find intermediate 
values 9n and r„ between r„ and r^^ such that 



where we recall ()3.12|) . This contradiction completes the proof of ()3.15|) . 
To complete the proof of the lemma we split the integral as 

r+{E,L) rrL rr+{E,L) 

= = / ...+ / ...=:/i + /2. 

_(E,L) j2E-2'^L{r) Jr.{E,L) Jr^ 



In the first term, we make a change of variables u = a/^E'l('") — ^l(^l) so that 

du 1 ^ 

dr 2u 



du_ i.^^(^)<Oon [r_(E,L),ri[. By 



J^e-^Utl) ^2{E-^>L{rL)-u^)du 

du /-^ ds 

< — / , = — / <oo 



by a further change of variables u= ^^/E^^WiJrL)s. The same type of esti- 
mate holds for the second part I2 of the integral under investigation, and the 
proof of the lemma is complete. □ 
We now show that the cut-off function hl^^ is square integrable and 
solves the equation {/o,/i1q,„} = (7lf7^ in the sense of distributions, more 
precisely: 

Lemma 3.3 For any mGN large, /il^^ G L^(M^), and for any spherically 
symmetric test function ip = il){r,w,L) G C"^([0,oo[x]R x [0,oo[), 



E<Eo J E<Eq 

Proof. We first prove that /ilf7,„ gL^(]R^). Since the integrand is even in v 
we can apply the change of variables ()3.14j) : 



In™ ht^dvdx 



ire' 



/ h\r,^2E-2^L{r),L) =dEdL, 

S^Jr-{E,L) ^/2E~2^L{r) 
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where 

Sm:=\{E,L) = {E,L){x,v)\{x,v)esnppfo, E<E^--, L>-\. (3.16) 
[ mm] 

Let {E,L) ^Sm- In the estimates below we write w{r) = v/2E-2^i(r) and 
r±{E,L) =r± for brevity. Then by the definition (j3.13p of h and the fact that 
l/\(j)',{E)\<lioT E<Eo, 



dr 

h'^{r,w{r),L)- 



w{r) 



< 



< 



r+ 



g{s,w{s),L) 



g{s,w{s),L) 



ds 



Wis] 



ds 

^ dr 



dr 
w{r) 



dr 



g {r,w{r),L) 



w{r) 
dr 



w[r ) 

^+2, ^ ~ dr 



wir] 



g (r,w(r),L)- 



wir] 



in the last two estimates we used the Cauchy-Schwarz inequality and 
Lemma 13.21 A further integration with respect to E and L and the change 
of variables {r,E,L) i— ^ {x,v) shows that ||/iln^||2 is bounded in terms of Cm 
and ||5'||2- 

Now let ip = ip{r,w,L) be a test function as specified in the lemma. Along 
characteristic curves of (j3.1ip . which as before we parameterize by r distin- 
guishing between w>0 and w <0, a simple computation shows that 

{E,iIj} = -signw^/2E - 2^^ (r ) ^ [V- (r, sign w - ir),L)]. 



By the change of variables used repeatedly above, using the abbreviation 
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dEdL 



w{r):= a/2£' — 2\l/i(r), and recalling the definition (|3.16p we find that 
{fo,ij}hln^dvdx= / ...+ / 

J{w>0} J{w<0} 

r r r+ d dr 

- (t>'o{E) w{r)-[^{rMr),L)]h{rMr),L)—-dEdL 

J JSm Jr- '^V) 

r r r+ d dr 

/ / / {-w{r))-[i,{r-w{r),L)]h{r-w{r),L)—- 

J JSm Jr. w{r} 

^;(E)£%.(r.-Mr).L) «';--';>;^> dr.EdL 

i)gln.^dvdx] 

notice that /i(r, ±w(r),L) = for r = r±{E,L), which together with the def- 
inition ()3.13p of h along the characteristics was essential in the integration 
by parts above, and also that g is even in v respectively w. The proof of 
Lemma f3. 31 is now complete. □ 

Regularization of hlfi^. 

The function hl^^ is not smooth, hence Lemma 11.11 cannot be applied to 
it, and therefore we smooth it. For fixed mGN the function hlQ^ is, as a 
function of r,w,L, supported in a set of the form 

Qm := [Ro,Ri] X [-Wo,Wo] X [Lo,Li] (3.17) 

with 0<Ro<Ri, Wq > 0, and < Lq < i^i; it will be important that its sup- 
port stays away both from r = and L = 0. Let C£C^(M'^) be even in all 
three variables, i.e., ^(2:1, 2:2, 2:3) = C(|2:i|, 1^21, ksl), C^O, J( = l, and define 
(n'-=n^C{i^') n^N. For n sufficiently large we define 



cx> poo 



hn{r,w,L) = / / / {hlny{r,w,L)(n{r-r,w-w,L-L)dLdwdr. 

Jo J~odJo 

Since is strictly inside the set {E<Eq} with a positive distance from its 
boundary, /i„ G C^(]0,oo[x]Rx]0,cxd[), and without loss of generality we can 
assume that supphnCQm^{E <Eo}. Since h is odd in w respectively v so 
is hn, and clearly, in L^flL^. The crucial step is to show that 

\im{fo,K} = gln^ (3.18) 
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in L^. To this end we fix {r,w,L), write E := E{r,w,L) and f := /o°°/^/o°° 
for brevity, and split tlie convolution integral into tfiree parts: 



{/o,/^n} = -0o(i5) 



wdr+(^^-U;,{r)^d^ 



Wdr + 



-,-Ul,{r)]d„ 



Cn{f — ,w — w , L — L) dL dw dr 



ihlnJif,w,L)<P',{E) 



= :-h + l2-h. 
By Lemma f3 .31 



wdr+\^-U'{f)\d, 



Cni'f' — r,w — w,L~ L)dLdwdr 



(nir — r,w — w,L — L)dLdwdr 



gln^{r,w,L)Cn{r-r,w-w,L-L)dLdwdr^glQ,^ in L^. 



We sliow tliat Ji and I3 tend to zero in L^. To do so we cliange tlie integration 
variables into r = n{r — f), w = n{w — w), L = n{L — L) so that f = r — r/n, 
w = r — w/n, L = r — L/n, dLdwdf = n~^dLdwdf, and drCni'f' — ^ — w , L — 
L) =n'^df({r,w,L) with analogous formulas for the other derivatives. Thus 



ihlnJ{r,w,L) 



1/n 



wdf + 



L 



3 f/o(r) 1 d,r, 



({r ,w ,L) dLdw dr , 



(hlaJ(f,w,L)4,'„(E) 



1/n 



1/n 



d,7, 



({r,w,L)dLdwd'r 



Now 



q^',{E)-^',{E) 
1/n 
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for n— s-cxo, and all these limits are uniform with respect to {r,w,L) EQr 
defined in (I3.17j) and with respect to (f,w,L) GsuppC- Hence 



h^-ihlnJir,w,L)<P'^iE) 
L rL 

= 0, 

h^ihlnJ{r,w,L)<p',{E) 
L 3Lf 



wdf+ { ^-U'oir) ] du 



({r,w,L)dLdwdr 



-wdf+ -^ + - 



4 u;;{ry 1 d,, 



({r,w,L)dLdwdr 







where both limits are in and the zeroes result from integration by parts 
with respect to f and w and an exact cancellation. The assertion ()3.18|1 is 
now established. 

Finally: The desired contradiction. 

The functions hn have all the properties required in Lemma 11.11 and we 
obtain the estimate 



D'ncifo)[{fo,hn}]>-- 



(j)o{E)-U;^{r)\hn\'^dvdx. 

{E<Eo} 1" 



By (I3.18P and the fact that /i„ and glnrr, are supported in a common compact 
set, V?7{/(,,/i„} — > Vf/c,if5^ in as n-^oo. Since ^U^ir) is bounded on [0,oo[ 
we obtain the estimate 

D'HcmalnJy-l [ <P'oiE)-Ul,ir)h'ln^dvdx. 



Since limm^oofl'ln™ =^7 in iv^nL^(M^), there exists mo sufficiently large such 
that s^lfirng 7^ and by Lemma IT!?! /iIq^^ ^ 0. For all m > mo, flmo C flm and 
/i^ln^^ < h^lQ^. Since 4>q{E) < and j.Uo{r) > we have for all m >mo, 



D'nc{fo)[glnJ> 



cP'^{E)-U'^{r) hHn^^ dvdx =: C > 0. 

{E<Eo} ^ 



Hence m — > oo leads to a contradiction to Lemma f3.H and the proof of The- 
orem |221 is complete. 
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4 Appendix: Proof of Lemma 11.11 



Let 



Uhix) :-- 



{fo,h}dv- 



dy 



\x-y\ 



be the potential induced by —{fo,h}, which is spherically symmetric. A short 
computation using the definition of the Lie-Poisson bracket shows that 

{fo,h}dv = Va:- j v(f)'Q{E)h{x,v)dv. 
By the spherical symmetry of Uh and h, 

J\x\<r 

1 



Kir) 



v(pQ{E)h{x,v)dv 



X 



^ J\x\=rJ r 
Therefore, by the Cauchy-Schwarz inequality, 



(f)Q{E)h{x,v)v ■ —dvduj{x) =4,71 j w(I)'q(E) h{x ,v) dv 



I \VMh\'dx<2n 

on 



- I w'^^o{E)dv 



- / (j)o{E)h^dv 



dx. 



Since 



an integration by parts with respect to w yields 



- / w'(l)'Q{E)dv- 



TX 



-'-oo 



0(0^ +^ + UQ{r)] dwdL = pQ{r) 



Hence 



D^'Hc{h)[{hM]>-\ j j <Po{E) {\{E,h}\'-Anpoh') dvdx. 



Since h is odd in v respectively w the function 



fi(r,w,L) := — h(r,w,L) 
rw 
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is smooth away from r = 0; in passing we notice that the functions /i„ to which 
we apphed Lemma 11.11 in the proof of Theorem 12.21 have support bounded 
away from r = 0, but this is not necessary for Lemma fl. II Since h = rwfi, 

{E,h} = rw{E,^} + ^{E,rw}, 

and hence 

\{E,h}\^ = {rwY\{E,fi}\'' + rw{E,rw}{E,fj.^} + ix''\{E,rw}\'' 

= {rwy\{E,iJ,}\'^ + {E,iJ,'^rw{E,rw}}-iJ,\w{E,{E,rw}}. 

The first term is as claimed in Lemma 11.11 The second term leads to 
{/o,//^rw{£',rw}} whose integral with respect to x and v vanishes; if we 
cut a small ball of radius e about x = from the x-integral then the surface 
integral appearing after the integration by parts with respect to x vanishes 
for e— i^O since rfi'^<C/r. By the Poisson equation, 

{E,{E,rw}} = -2wUo-w{rUQy = -rw (^Anpo + ^U;^^ . 

Hence the third term above becomes Airpoh'^ + h'^^UQ, and Lemma |l.ll is 
proven. □ 
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